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Abstract 

It is shown that the continuous g-Hermite polynomials for q a root of unity have 
simple transformation properties with respect to the classical Fourier transform. This 
result is then used to construct g-extcnded eigenvectors of the finite Fourier transform 
in terms of these polynomials. 

2000MSC: 33D45, 42A38 

1 Introduction 

The finite Fourier transform [26\ Ch. 7] (also called discrete Fourier transform) is denned 
as the Fourier transform associated with the finite abelian group Z^r := Z/(iVZ) of integers 
modulo N \25\ 119]. just as the classical integral Fourier transform is the Fourier transform 
associated with R. In concrete terms, it is a linear transformation $W of the space of 
functions on Z with period N defined by 

($W/)(r) : = _ (-£ raj f(s), r G Z. (1.1) 

Equivalently, if we identify the iV-periodic function / on Z with the vector 

(/(0),/(l),...,/(iV-l)) 
m then $^ is a unitary operator on C N with matrix elements 

*W = exp f ^rs J , 0<r,a<N-l. (1.2) 



We are interested in a suitable basis of N eigenvectors f^ (n = 0, 1, ...,N — 1) of 
which thus should satisfy the eigenvalue equation 

2V-1 

E^ ) fk N) (s) = X n fi N \r) (1.3) 

s=0 
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for suitable eigenvalues A n . Since the fourth power of $>( N ) is the identity operator (or 
matrix), the A n 's can only be equal to ±1 or ±i. 

The finite Fourier transform has deep roots in classical pure mathematics and it is also 
extremely useful in applications. See [12], [27] for mathematical and historical details and 
|13j . [14] for the relation with Gauss sums. 

Mehta studied in [22] the eigenvalue problem (|1 .3|) and found analytically a set of 
eigenvectors Fn of the finite Fourier transform $W of the form 

oo 

Fi N) (r):= £ e~^ kN ^ 2 H n (^ (kN + rj), r = 0, 1, . . . , N - 1, (1.4) 

k=— oo 

where H n (x) is the Hermite polynomial of degree n in x. These eigenvectors Fn ^ corre- 
spond to the eigenvalues A n = i n , that is, 

N-l 

£ *W =i»Ff)(r), r = 0, . . . , iV - 1. (1.5) 

They can be considered as a discrete analogue of the well-known continuous case where the 
Hermite functions e~ x I 2 H n {x) are constant multiples of their own Fourier transforms: 

/ j*v-*V* H n (x) dx = i" e-y 2 ' 2 H n {y). (1.6) 

It was conjectured by Mehta [22] that the (n = 0, 1, ... ,N — 1 if TV is odd and 

n = 0, 1, . . . , N — 2, N if A^ is even) are linearly independent, but this problem has not 
been any further considered until now. As shown by Ruzzi [23], these systems are in 
general not orthogonal. 

In later work it was shown that many q-extensions of the classical orthogonal poly- 
nomials satisfy simple transformation properties under the Fourier transform (see [7J and 
references therein). Thus it was natural to repeat Metha's construction of eigenfunctions 
of the finite Fourier transform in that context. This was done in [S], [8] and [11] . In 
particular, it was shown there that the finite Fourier transform provides a link between 
continuous g-Hermite and g _1 -Hermite polynomials of Rogers, as well as between families 
of Rogers-Szego and Stieltjes-Wigert polynomials. It turned out that the same form of 
connection exists also between discrete g-Hermite polynomials of types I and II, see [3]. 

Zhedanov [28] considered continuous (/-Hermite polynomials for q a root of unity and he 
obtained a discrete orthogonality on finitely many points with complex weights for a finite 
system of such polynomials. In the present paper we consider the transformation properties 
of these polynomials under the integral and finite Fourier transforms. Analytically, this 
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turns out to be a straightforward extension of the earlier results for the continuous q- 
Hermite polynomials with < q < 1. However, the resulting formulas, a little different 
from the case that q is real, are interesting enough to be displayed. 

Another feature of the present paper, compared with [SJ, [8] and [11] ; is that we 
emphasize a more conceptual approach by using Theorem 14.11 due to Dahlquist [16] and 
Matveev |21j and the (trivial) Lemma |5.1[ rather than repeating a technical argument in 
each special situation. 

The contents of the paper are as follows. In section 2 we recall some properties of the 
continuous g-Hermite polynomials for general complex q, and in particular for q a root 
of unity. In section 3 we consider the behaviour of these polynomials times a Gaussian 
under the integral Fourier transform. This result, together with the Dahlquist-Matveev 
Theorem l4.lt then gives a construction of functions behaving nicely under the finite Fourier 
transform. In section 5, using Lemma [5. II we obtain from these functions and their Fourier 
images eigenfunctions of the finite Fourier transform. Finally, section 6 concludes the paper 
with a brief discussion of some further research directions of interest. 

Throughout our exposition we employ standard notations of the theory of special 
functions (see, for example, [2] and [TT]). 



2 g-Hermite polynomials for q a root of unity 



The continuous g-Hermite polynomials of Rogers (see [23], [I], [6]), denoted by H n (x\ q), 
can be generated for any q £ C by the three-term recurrence relation 



2x H n (x\ q) = H n+1 (x\ q) + (1 - q n ) H n _ x (x\ q) 



(2.1) 



with initial condition Hq (x\ q) = 1. Their explicit form as a finite Fourier series in 
(x = cos 6) is given by 

n i- -i 

H n { C os9\q)=Y J " 

k=0 



k 



j.(n-2k)9 



(2.2) 



where the symbol [^] stands for the q-binomial coefficient 



{q;q)n 



(q;q)k (q;q) n -k 



n 

n — k 



(2.3) 



Here (a;q) n is the g-shifted factorial, see [TTJ (1.2.15)]. The right-hand sides of (j2.3|) and 
(|2.2p are well-defined for all q £ C because the (/-binomial coefficients are polynomials in q. 
Since sin# = cos(|7r — 6) we can rewrite (|2.2p as 



H n (sin ( 



E 



k=0 L j 1 



k i(2fc-n) 



■lfe 



(2.4) 
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The polynomials H n {x\ q) are orthogonal polynomials for < q < 1. For q > 1 they 
are orthogonal polynomials in ix : the g -1 -Hermite polynomials (see [5]) denoted by 

h n {x\q):=r n H n {\x\q- 1 ). (2.5) 

Another case of orthogonality, but not with positive weights, was considered for q a 
root of unity, see [28J. For M a positive integer put 

^• M :=exp(27rij/M), j G {1, 2, . . . , M - 1}. (2.6) 

For such q = qj t M (|2.ip and (|2.2p remain valid. In particular, for q = q^u with j and M 
co-prime and for n = M, the only non-vanishing terms in (|2.2p occur for k = and M. 
Hence, for j and M co-prime, we have 

H M (cos9\ q jiM ) = 2 cos M(9 =: 2T M (cos0) , (2.7) 

where Tm(x) is a Chebyshev polynomial of the first kind (see, for example, j2j Remark 
2.5.3]). 

As pointed out in [28], the polynomials H n (x\ c/^m) (n = 0, 1, ... ,M — 1) satisfy a 
discrete orthogonality with possibly complex weights on the M zeros of Tm{x) if j and 
M are co-prime. It is for functions suitably defined in terms of these polynomials that we 
will discuss their integral and finite Fourier transforms. 

Finally observe that, by induction with respect to m and n, we derive from (|2.ip and 
p. 71) that, for j and M co-prime, 

H mM +n(x\ qj,u) = (2T M (x)) m H n (x\ q j)M ), n = 0, 1, . . . , M - 1, m = 0, 1,... . 

(2.8) 



3 Integral Fourier transform 



There are g-extensions of the eigenfunction result (|1.6j) for the integral Fourier transform. 
These interrelate certain ^-polynomial families (see [7] and references therein). For the 
continuous g-Hermite polynomials we obtain: 

Lemma 3.1. The Fourier transform of the functions e~ x2 ^ 2 H n (sm(Xx)\ q) is given by: 

I POO 

J_ / e i X y- x 2 /2 Hn (sm(\x)\q)dx = i n e- n2x2 / 2 e-y 2 / 2 

V27T J-oo 



E 

fc=0 



q^e^ 1 



)f^{ ^ j 
(-l) fc e" (2fc - n)Ay , A,g G C, q ^ 0,1. (3.1) 
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In particular, if q = e 2x2 then 



-i roc 

-== / e ix y- x2 / 2 H n (sm(\x)\q)dx = q n2 / 4 e-y 2 / 2 H n (sm(i\y)\q~ 1 ). (3.2) 

V27T J-oo 

Proof Substitute (|2.4|) on the left-hand side, take termwise Fourier transforms (which 
turns down to the Fourier transform of the Gaussian), and use that 



k(k—n) 



□ 



If q = e ~ 2x2 then we see by ([23]) that ([33]) simplifies to (l33]h 

We consider two cases of equation ()3.2[) . First let < q < 1 and define k by 

g = exp(-2K 2 ), 0<k<oo. (3.3) 
Then, by ([23]), we obtain (see [TO]. [7]^: 

Proposition 3.2. The Fourier transform of the functions e _x ' 2//2 i? n (sin kx\ q) is given by: 
i r°° 

e ™y-* 2 /2 Hn (smKx\q) dx = i n q n2/4 e- y2/2 h n (smhKy\q). (3.4) 

V27T J -oo 

Second, let q := q^M as in (|2.6[) and put A := ctj^M, where 



-4 e m//4 , hence e 2Qj ' M = 7 - m and \ctj m = ctj M ■ (3-5) 



Proposition 3.3. The Fourier transform of the functions e x2 / 2 H n (sm(aj^M x )\Qj.Ai) is 
given by: 

—7= [ e lxy ~ x2/2 H n {s\n(a j)M x)\ q jiM ) dx = o n J 4 e~ y2/2 H n ( sin(a]Jiy)\ qj^) . 

(3.6) 

The Fourier inversion formula of (|3.6|) is just the result of taking complex conjugates 
on both sides of (|3.6p . We will mostly work with (|3.6p for n = 0,l,...,M — 1, but this 
formula remains valid for all nonnegative integer values of n. In particular, for n = mM 
(m = 0, 1, . . .), formula (|3.6|) takes by substitution of (|2.8p the form 



1 



cos 



(Mtt/2- ^jMe- if / 4 x) e ^ 3:2/2 dx 

= i m2j ' M e" y2/2 cos m (Mtt/2 - e^ 4 y) . (3.7) 



Depending on the value of M (mod 4) this may be further simplified. 
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Remark 3.4. Formulas (|3.4p and (|3.6|) can be considered as g-analogues of formula (|1.6p . 
in the sense that (jl.6p can be obtained as the limit for q f 1 of (|3.4p and as the limit for 
M — > oo of (|3.6p with j fixed. Indeed, from |18^ (5.26.1)] we have 

lim - q) y U H n (x^(l -q)\qj = H n {x). (3.8) 



Since k ~ ((1 — q)/2) 2 as g f 1, it follows from (|3.8p that 

lim K~ n -ff n ( sin(«;x)| q) = H n (x), 
<?Ti 

hm (WT n il n ( sin(i/cj/)| g _1 ) = H n (y). 
Hence, in view of (|2.5p . equation (|3.4p with both sides multiplied by n~ n tends to (jl.6p as 

ffTi. 

As for (|3.6p with j fixed, we have ol^m ~ ((1 — Qj,M)/^) 2 as M — > 00. Hence it follows 
from (OLED that 



lim a ■ n M H n (sin(atj t M x)\ q jiM ) = H n (x), 
lim F n (sin(5^ y)\ qAj) = H n (y). 

Then (|3.6p with j fixed and with both sides multiplied by aj^, tends to (jl.6p as g | 1. 

4 Finite Fourier transform 

It was observed by Dahlquist [16, Theorem 1] and Matveev |21l Theorem 8.1] that Mehta's 
result (jl.5p is a special case of the following more general relationship between integral 
Fourier transform and finite Fourier transform, which can be obtained as an immediate 
consequence of the Poisson summation formula. This result will hold for / in a wide 
class of functions on R, but for convenience we only formulate it for / £ S, the space of 
Schwartz functions on R (see [26, Ch. 5, §1.3]). 

Theorem 4.1. Define a linear map A4^ N ' from S to the space of N -periodic functions 
on Z by 

( M Wf)(r):=J2f(\/¥(kN + r)), feS,rGZ. (4.1) 

kez 

If f,g 6 <S are related by the integral Fourier transform 

g(y) = -L f°° e i ^f(x)dx, (4.2) 
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an 



d if F := M.( N >f, G := A4^ N 'g then F and G are related by the finite Fourier transform 



N-l 



G(r) 



s=0 



2vri 



exp [ rs | F(s). 



(4.3) 



In particular, in the case when g{x) = A/(x), A = ±l,±i, one has G(r) = XF(r); so 
Mehta's eigenvectors (jl.4p are a particular case of the general statement. 

Let us now apply Theorem 14.11 to the case where / and g are implied by (|3.2p . i.e., 

f{x) := e -* 2 / 2 H n { S m(\x)\q), g(y) := e^ 2 F n (sin(iAy)| g" 1 ), 

where A, q € C, q ^ 0, 1, and q = e~ 2x2 . Put 

oo 

fi N \r\q):=(M^f)(r)= £ e~f (sin (a^ + r)) | q) , 

k=—oo 

9 { n N) (r\q) := q -" 2 /*(MWg)(r) 

oo 

= e'ff^+^H^sm^X^ikN + r^q- 1 ). 

k=— oo 

Then we obtain by Theorem 14.11 

Proposition 4.2. The finite Fourier transform of the functions fn N ^ is given by: 



(4.4) 



(4.5) 



2V-1 



s=0 



*™ 1 ?))(r) = £ exp ( ^r, ) f^(s\ q) = q^g^(r\ q). (4.6) 



N 



Just as we did in §3 for (|3.2|) . we consider two special cases of (|4.6p . First, let < q < 1 
and take A = k as in (|3.3|) . Then (|4.6I) holds (see [9] and [TT]) with 



oo 

fl^°(r|9)=i B E e-*< wv ^X(BiBli(K V /^(fcJV + 



(4.7) 



fc=— oo 



Second, consider (|4.6|) with g = q^j^ and A = Oj,m as in (|3.5|) . Then (|4.6|) holds with 



oo 

/r ) H9i,M)= E e-^( fc ^) 2 ^(sin(a,, A/ yf (fciV + 

fc=— oo 

oo 

fl^(r|<fc*r) = E e-#^) 2 ^ t (sin(a^7yf (fciV + 



fe=— oo 



*3,M 



(4.8) 



(4.9) 



7 



By similar reasoning as in Remark 13.41 we see that (jl.4|) is a limit of (|4.6p in these 
two cases (as q f 1 in the first case and as M —* oo in the second case). These limits are 
formal because we have to take termwise limits for fn\ s \ l) an d gn{ r \ <?)• Just as for 
the functions Fn N ^ given by (jl.4p . we do not know if the functions fn defined by (|4.4p 
are linearly independent. 

Let e = ±1. From (|4.ip we see that, if /(— r) = e/(r) for all r, then also f)(— r) = 

e(.M (Ar) /)(r) for all r. Also note from d^TJ) that F n (-x| <?) = (-l) n # n (a?| </)• Hence, by 
(lOP and (03}, 

/^(-rk) = (-ir/^(r|g), ^(-rjg) = (-l) n gW(r\q). 

Thus consider (|4.6p with /i^' and #n as in (|4.8p . (|4.9p . Then take complex conjugates 
on both sides and use the simple facts just mentioned above. Then we obtain: 

* (A ° (A N) ( • I QJM)) (r) = %m & N) (r I ^,m) , (4.10) 

* (JV) (/l*°(.|^))(r) = (-ir^^/^Crl^M). (4.11) 

Remark 4.3. It is tempting to consider the case q = q* m an d ^ = a j,M of (|4.6p with 
M = N, even if one loses then the possibility to take the limit to (|1.4p for M — > oo. 
One might hope to arrive at some transform acting on the polynomials H n {x\ CRjv) { n 

0. 1 \ 1; j and -/V co-prime) but only involving their values at the zeros of T/v(x), 

1. e., at the points involved in Zhedanov's [28J orthogonality relations for these polynomials. 
At the moment we have no idea how to proceed here. 

5 Finite Fourier g-extended eigenvectors 

We now show that relations (|4.10p , (|4,lip enable us to construct ^-extensions of Mehta's 
eigenvectors (|1.4p of the finite Fourier transform when the deformation parameter q is a 
root of unity. For this we need the following trivial observation. 

Lemma 5.1. Let V be a complex linear space. Let f,g £ V and a,b £ C such that 
= a 2 g and <&g = b 2 f. Then 

^(bf±ag) = ±ab(bf±ag). 

Combination of this lemma with (I4,10|) and (|4,lip shows that the functions 

■ n -n 2 /8 f (N)f I \ i n 2 /8 AN), , 7 

1 Qj,N fn '(■\Qj,M) ±( lj,M fn ( • I Qj,M) 
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are eigenfunctions of with eigenvalues ±i n . So put 



= «(e tatl/8 rf /8 /^(r|( & ,M)) 
= 9( e iW8 ^M /8 e ) (rk J ,M)) 



Then 



g ij M))(r)=i"FW(r|^ )M ) ) 
g il M))(r) = -i n GW(r|g i>M ). 



(5.1) 
(5.2) 



6 Concluding comments and outlook 

We have demonstrated that the continuous (/-Hermite polynomials for q a root of unity are 
interrelated by the classical Fourier transform, see (|3.6|) . Then the technique of construct- 
ing eigenvectors of the finite Fourier transform, pioneered by Mehta [22] and formulated 
in a more systematic way by Dahlquist [16] and Matveev [21], has been employed in order 
to construct (/-extended eigenvectors of the finite Fourier transform. 

Quite naturally, it would be of considerable interest to find out whether there are other 
families of (/-polynomials for q a root of unity, which also possess such simple transforma- 
tions properties with respect to the Fourier transform and, consequently, give rise to other 
(/-variants of Metah's eigenfunctions (jl.4p of the finite Fourier transform. The point is 
that the continuous (/-Hermite polynomials, considered in the present paper, belong to the 
lowest level in the Askey hierarchy of basic hyper geometric polynomials [18]. Therefore it 
will be natural to attempt to apply the same technique to the study of other (/-families, 
which depend on some additional parameters (and therefore occupy the higher levels in 
the Askey (/-scheme). 

Finally, another direction for further study is connected with (/-extensions of the har- 
monic oscillator in quantum mechanics [3], [20], [15]. We remind the reader that for 
proving the fundamental formula (jl.5p for Mehta's eigenvectors F^ n ' of the finite Fourier 
transform operator (II. 2D it is vital to use the simple transformation property (|1.6p of the 
Hermite functions H n (x) exp(— x 2 /2) with respect to the Fourier transform. Moreover, 
these eigenvectors are actually built in terms of these Hermite functions taken at the 



infinite set of discrete points xf> := J ^ (kN + r) , < r < N - 1 , k G Z (cf. flgD ). 
In other words, Mehta's technique of proving (jl.5p is based on introducing a discrete ana- 
logue of the quantum harmonic oscillator. It seems that (/-extended eigenvectors of the 
finite Fourier transform, constructed in the foregoing sections, can be similarly viewed as 
discrete analogues of the (/-harmonic oscillator of Macfarlane and Biedenharn for q a root 
of unity. 
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We plan to continue our studies in both of these directions. But a focus of our at- 
tention will be on constructing an explicit form of a finite-difference equation for suitable 
eigenvectors fjl .3|) of the finite Fourier transform operator $( Ar ). From the group-theoretic 
point of view this amounts to finding an adequate finite representation, associated with the 
Heisenberg-Weyl group [25l [19] . Purely analytically this is reduced to the construction of 
an explicit finite-difference operator L which commutes with so that the eigenspaces 

of $W can be decomposed by spectral decomposition of L. 
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